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ABSTRACT 


The objective of this study is to model and simulate using an Analog 
computer the performance of a relatively complex engineering system. The 
results of the simulation are validated by a comparison of the simulation 
with the measured performance of the prototype system. 


The prototype system consists of a single shaft gas turbine engine 
driving a three-phase alternator. The alternator load is a variable speed 
three phase wound rotor induction motor which drives a ship bow thruster 
propeller. In order to provide sufficient data to validate the simulation, 
the alternator speed, current and terminal voltage, KW load, and induction 
motor speed are recorded continuously as the motor is switched through five 
speeds. The recorded transients in these variables are then used to validate 
the model and simulation. 


The system is segregated into five sections and analytical models are 
determined for each component. The gas turbine is modeled using a transfer 
function derived from the engine performance curves. The speed governing 
System is modeled using transfer functions for the governor components zs 
Supplied by Woodward Governor Company. The alternator model is a result of 
the work of Professor Paynter and is taken from notes from his subject 
"Dynamics and Control of Rotating Machine Systems." The induction motor 
model is derived from a simplified equivalent circuit. The propeller is 
modeled as a torque proportional to the square of the speed. The simula- 
tion is simplified considerably by using a single-phase model to represent 
the electrical machinery. 


The correlation between the simulation and the prototype performance 
is excellent. The models used are relatively simple and yet can provide 
close estimates of the system performance and can demonstrate how perform- 
ance varies as the system parameters are changed. 
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Thesis Reader: Warren C. Dietz 
fille: Associate Professor of Marine Engineering 





2. 


ACKNOWLEDGEMENTS 


I would like to first thank my wife Diana for the patience she has 
shown during the last two years of school work, and especially during the 
last term. 

During this study valuable assistance has been received from many 
individuals and in various ways. The most important has been the advice 
and guidance received from Professor Paynter, both during this study and 
in previous course work. 

I am particularly grateful to Mr. Henry Egan of the Boston Naval 
Shipyard for providing me with the equipment required for monitoring the 
performance of the prototype systen. 

The cooperation of the United States Coast Guard First District 
Naval Engineering Branch and the engineers of the USCGC BOUTWELL and the 
USCGC CHASE was indispensible in making the required recordings of the 
system operation. In particular, the willingness of Chief Electrician 
Duke to allow an outsider to connect wires and recorders to various Sec- 
tions of the BOUTWELL'S electrical equipment was most encouraging. With- 
out his willing help, it would have been most difficult to obtain system 
data. 

Also, I would like to thank Richard S. Sidell for his assistance 
over the past four terms with the EAI-680 Analog computer. 

Finally, I would like to thank Mrs. Jack "Sam" Price for her excel- 


lent and efficient work in typing this paper. 





pa 一 em zx Am ze mm om‏ — = — + ہے 
v ÄR - —‏ - عم ہے 
su a. ami a n‏ = - س سد سد 4 < 


— Á —r سےا‎ - —  — — 


m 
o 
(n 
ct 
O 
rh 
۲ 
ja 
pa 
ke 
LU) 
it 
> 
fu 
rt 
Lie 
0 
= 
0۷۸ 
y 
N 
(Y 
ej 
b; 
y 
1 一 * 
(b 
Uu 


D. 


y. 


Refe 


instrumentation and Recoräin 
System Models 

4.1 Gas Turbine 

4.2 Governor Section 

1.3 Alternator 

BD dnduction Motor 

۲۰٠٦ 0876661167 

4.6 Bearing Resistances 
Simulation 

5.1 Parameter Estimates 

5.2 Analog Setup 

553” Sealiag Analog Simulation 
Results 


Cerelusions and Recommendations 


Bences 


vu 
íM 
09 


(D 


N) 


(1) 


pa 
(a) LO) o (J 


ho 
(1) 


{> (1) (1) 
pa N3 e» 


$ 
rn 


Dk 
O) 


ER 
Í 


U1 
O 


O) 
~ 


00 
I 


OO 
10 





Ens 


JAR 


T2. 


13. 


ce 


LIST OF ILLUSTRATIONS AND TABLES 


Illustrations 
Prototype System Bond Graph 
Cutaway View Saturn Gas Turbine Engine 
Reduction Gear Drive Arrangement 
Saturn Engine Performance Curves 
Governing System Schematic 
Schottel Bow Thruster Arrangement 
Propeller with Nozzle 
Instrumentation Bond Graph 
Analog Computing Diagram 
Prototype Measurements and Simulation Results 
Induction Motor Performance Curves with Simulation Results 
Engine and Governing System Block Diagram 


Engine and Governing System Block Diagram 


Tables 
Amplifier Assignment Sheet 
Potentiometer Assignment Sheet 
Correlation of Alternator Speed 
Correlation of Induction Motor Performance 


Comparison of Estimated Parameters with Values Used in Simulation 


Page 
16 
17 
18 
19 
20 
21 
22 
27 
57 
69 
83 
31 


28 


64 
66 
78 
19 


80 





ALT 


a9 


は 3 


F CV 


GCB 


-6- 


7  + 5 


Alternator 

Governor constant 

Fuel control actuator buffer spring stiffness 
mueci=cContrel actuator constant 

Fuel control actuator needle valve setting 
Propeller water velocity angle 

Propeller torque coefficient 

Current transformer: 

Propeller diameter 

Voltage 

Governor error signal 

Output of governor amplifier section 

Fuel control actuator 

Fuel control valve 

Governor control box 

Governor stabilizing section damping ratio 
Alternator line current 

Alternator field current 

Inertia of all components driven by gas turbine 
Inertia of motor 

Inertia of propeller 

Induction motor 

Gain of governor amplifier section 


Fuel control actuator constant 





max - 
min 


A | 


Engine torque constant 

Alternator inductance 

Motor inductance 

Speed in RPM 

Deviation of alternator speed from base speed 

Base speed of alternator and motor 

Air flow rate 

Fuel flow into engine manifold 

Fuel flow into fuel control valve 

Fuel flow into engine manifold 

Ambient pressure 

Engine exhaust pressure 

Potential transformer 

Bearing resistance of engine and alternator 

Bearing resistance of induction motor and propeller 
Motor rotor resistance 

Added rotor resistance 

Dynamic operator representing fuel control valve action 
Water density 


Laplacian operator 


Angle of fuel control actuator shaft 


Maximum travel of fuel control actuator shaft 
Engine mechanical time constant 
Torque 


Normalized torque (T/Tyer) 





Alternator torque 

Engine torque 

Motor torque 

Engine inlet temperature 

Engine exhaust temperature 

Bond Graph symbol for engine speed reducer 
Bond Graph symbol for motor speed reducer 
Axial water velocity through propeller 
Tangential velocity of propeller blade at maximum diameter 
Speed in radians per second 

Motor synchronous speed 

Electrical reactance 

Equivalent reactance of motor 


Synchronous reactance of alternator 





+92 


1. INTRODUCTION 


Simulation has become one of the most important tools available 
to engineers in the study and design of complex physical systems. The 
ability to predict a system's response to various inputs is useful in 
many ways. Actions which can have disastrous effects on a system can 
be avoided if accurate predictions of the system response are avail- 
able. Simulation can also aid in the realization of high levels of 
system performance by aiding in control system design. Studies of the 
effects of varying system parameters can be performed with a simulation 
so that an optimum or near optimum prototype can be built without a 
costly trial and error process. 

The first stage in simulating a complex physical system is to i- 
dentify the individual components. An attempt is then made to formulate 
a model for each component of the system so that the interconnected mod- 
els will predict the performance of the system within a specified level 
of accuracy. Simulation is in effect a four-stage process. In the 
first stage the system is reduced to individual components and "black 
boxes" are used to represent these components. In the second stage, 
the components are interconnected so as to reflect the structure of the 
prototype system. The third stage usually requires the use of a com- 
puting machine to determine the time history of the variables as the 
components interact with each other. 

This simulation procedure can be carried out almost completely us- 
ing the concept of Bond Graphs (1), (2). Although bond graph techniques 


can be used to develop a set of state equations representing the entire 
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system, for this simulation bond graphs will be used primarily to dis- 
play the structure of the system. The Bond Graph method enables one 
to visualize a system as a set of multi-port devices connected by power 
bonds. As a result of this partitioning of the system the energetic 
interactions across the boundaries are clearly shown. For the case of 
the prototype system under study the basic bond graph representing the 


main power flow is as follows: 


Gas Near EN Tnduüuction LN Propeller 


Turbine Motor 


Each component is a multiport which interacts with the rest of the 
system through the power bonds as shown. Associated with each bond 
are the power variables (T, Uu or e, i) which determine the power flow 
between the multiports. 

The models required for the multiports shown above can be deter- 
mined in various ways. However, regardless of how the models are de- 
termined, they should satisfy the following conditions: 

1) The time history of the power variables at each port when connected 
to the system should in the simulation match that of the corresponding 
prototype parts over the ranges of variables and frequencies being con- 
sidered. 

2) Each model should be as simple as possible consistent with the ac- 
curacy required for the simulation. Unnecessary complexity will in- 
crease the cost of simulation and will tend to obscure the more impor- 
tant points. More important however is the need to adapt the simulation 


to available computer facilities. Frequently, as in this case, the 





E 
simulation must be simplified because of limitations of the available 
computing equipment. 

There is no unique method of modeling system components. However, 
any method can be classified as being either analytical or empirical. 
Analytic methods rely on use of the basic engineering disciplines and 
in effect derive a model from first principles. Although in principle 
analytic methods can be used to model any device, these methods become 
extremely difficult to use for very complex components. Empirical 
methods are in effect a transfer function approach to modeling. Given 
a component and information about its performance under specified con- 
ditions, it is usually possible to infer a relationship between the in- 
put and output variables. For multiple input-multiple output systems 
this technique may become difficult to use. Also, this method gives no 
information about the internal operation of the component. However, for 
studies where only interactions with the rest of a system are considered 
this method is perfectly adequate. 

In this study the models are determined using both methods. The 
alternator model is based on an analytic derivation whereas the model 
for the gas turbine is derived from performance curves. The models used 
for this simulation are discussed in detail in the following sections. 

The question most frequently asked during this study was "Why simu- 
late a system which has already been built?" The final stage in any sim- 
ulation is validation of the results. Although in most cases the simu- 
lation must be validated without prior knowledge of the prototype, models 
which have been proven with one prototype can be used with great confi- 


dence for a large family of similar systems. In effect, the best 
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documentation of a system model is by comparison with observed perform- 
ance of a prototype. A large part of this study deals with collecting 
sufficient operating data from the prototype system in order to validate 
the models used. 

The prototype was selected because it provided the opportunity to 
model an interconnected system which included mechanical and electrical 
machinery with feed-back control on the energy source. Equally im- 
portant was the ability to instrument the system in order to record 
sufficient operating data to validate the model. The prototype selected 
was the gas turbine emergency generator set and bow propulsion system on 
the USCGC BOUTWELL which is located in Boston, Massachusetts. The pur- 
pose of this study is to simulate the effects of the transients in the 


system caused when the bow thruster motor is started. 
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S SEP CEI TITO OPI PPOTOTYPE SYSTEM 


The prototype system is shown in detail in the bond graph shown 
in Figure 1. The prototype system is defined from an arbitrary choice 
of equipment installed on the ship. The emergency generator can be 
used to supply electrical power to any load on the ship as shown in the 
sketch below. The bow thruster is a trainable propeller which is used 
to provide a directional thrust at the bow when docking the ship. The 
unit is retracted into the ship when not in use. The electrical system 
can be interconnected so that the bow thruster is the only load on the 
emergency generator. The prime mover is a single shaft Solar Saturn 
Gas Turbine. For standard atmospheric conditions the engine is rated 
at 1200 HP at 22,300 RPM. The output is reduced through a two stage 
planetary gear box to 1200 RPM. An engine cutaway view, arrangement of 
reducing gear and engine performance curves are given in Figures 2, 3, 
and 4. 

The engine speed is controlled by a Woodward EG series governor. 
The governor system is shown in the schematic on page 20. The two main 
components are the EGA Control Box and the EGB-2C Actuator. The gover- 
nor can control the engine in either isochronous or droop mode. The 
EE EE can also control the speed using the conventional mechanical 
flyweight assembly. For this study SR system will be controlled 
electrically and in the isochronous mode. The cperation can best be 
described by the transfer functions for the components. The transfer 
functions supplied by Woodward Governor Company are of the following 


form: 
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These signals are shown in the governor schematic, Figure 5. 

The error signal is a combination of a voltage proportional to the 
error in speed and a voltage output from the load sensing network which 
provides lead compensation. The governor response can vary over a wide 
range depending on the values of the following parameters: 

(1) Feedforward gain of the load sensing network (0.00405) 

(2) Gain of the amplifier section (K) 

(3) Damping ratio of the stabilizer section (y) 

(4) Buffer spring stiffness (a) 

(5) Needle valve setting (a3) 


The alternator is a General Electric machine with the following 


characteristics: 
Rated KVA 625 Speed 1200 RPM 
Rated KW 500 Winding 3 Phase 
Rated P.F. 0.8 Connection Y 
Rated Volts 450 ES SONES 
Rated. Amps 801 Poles 6 


Voltage regulation is provided by a static excitation system. 
The alternator is connected to the ship's electrical system in the 


following way. 
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For this study the system is set up so that the emergency generator is 
supplying only the Bow Thruster load. 

The Bow Thruster Unit is shown in Figures 6 and 7. It is manufactured 
by Schottel of America, Inc. The driving eer is a General Electric three- 


phase wound rotor induction motor with the following ratings: 


Rated Volts HHO Rated Amps 440 
Speed 1180 RPM Torque 1550 1b-ft 
HP 350 Frequency 60 Hz 


"The motor performance curves are shown in Figure ll. The motor is started 
with full line voltage applied and additional resistance in the rotor cir- 
cuit. Five propeller speeds are provided by shunting the external re- 
Sistors in four steps. The shunting of these resistors is controlled by 
time delay relays. The motor is connected to the propeller through a 
vertical shaft and bevel gears. 

The propeller is enclosed in a nozzle as shown on page 22. The 
diameter is four feet, pitch 28.819 inches and the reduction ratio from 


the motor speed 2.079:1. 
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Bow Propeller With Nozzle 
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3. INSTRUMENTATION AND RECORDING 


Verification of the system model and simulation by comparison with 
the prototype performance requires records of the system variables during 
operation. The optimum choice for recorded variables are the state varia- 
bles of the system model. In this case the choice of observed variables 
was limited to measurements which could be made without major modifica- 
tions to the prototype and which did not require excessive time or cost. 
As a result, most of the observed variables were chosen because the proto- 
type was instrumented to monitor these variables. However, these con- 
straints did not prevent obtaining a set of observables sufficient for 
model verification. 

The complete instrumentation setup is shown in the bond graph in 
Figure 8. The main power flow through the system is shown with heavy 
lines. The boxed items show all sensors, signal conditioners and record- 
ers which were added to the prototype instrumentation. Each box is la- 
beled with the variable recorded. All bonds to the boxes are active, 1. 
e.„ zero power flow into the recorders is assumed. Details of each meas- 
urement is given below. 

Results were recorded on a Clevite Mark 220 Brush Recorder. This 
instrument is very versatile and provided satisfactory records of all 
variables measured. Input signals ranged from five millivolts to fifty 
volts and in frequency from D.C. to 60 Hz. For full scale deflections 
the frequency response is flat up to 40 Hz. The frequency response can 
be increased to 100 Hz if the meter movement is limited to twenty percent 
of the full scale. 


ALTERNATOR LOAD The EGA Control Box Load Section senses the 





= Olt 
alternator terminal voltage and current through potential and current 
transformers. A correction is made for Br factor so that the output 
of the load section is a signal with a D.C. level proportional to the KW 
load. At full load the D.C. level is approximately nine volts. However, 
this value will vary depending on the values of the resistors used with 
the current transformers, the levels of voltages and currents delivered 
by the instrument transformers and the gain of the load sensing Section. 
The load signal contains harmonics of the line frequency in addi- 
tion to the D.C. level. Since the amplitude of these harmonics was ap- 
proximately twenty percent of the D.C. level an RC low pass filter with a 
time constant of 0.068 seconds was used to filter the harmonics. 
>> FIELD CURRENT The full load field current is approximately seventy 
amps at seventy volts. A shunt resistor was inserted in the field circuit 
at the brush rigging. The shunt resistor voltage drop is fifty millivolts 
at m rive amps. | 
ALTERNATOR SPEED The installed tachometers and frequency meter are 
used primarily for steady state observations. Large errors in measuring 
transients can be expected when the time interval of the transient ap- 
proaches the time constants of the meter movements. As a result, the 
meters cannot be used to accurately record transients resulting from step 
changes in load. 
The speed measurement was made using a Clevite Corporation Frequency 
Converter. The converter is calibrated to give a D.C. signal proportional 
to the change in frequency. A deviation of + 5 Hz from the base frequency 


gives a + 2.5 volt output. The rise time of the output signal is 0.060 





ےت 
seconds. This response together with the response of the recorder gives‏ 
an accurate measure of the speed transient.‏ 

The signal input to the converter is at 110 volts in the range 55-65 
Hz. The prototype system was set up with the alternator supplying the 
emergency switchboard so that the bow thruster was the only load. In 
order to provide a 110-volt signal at the alternator frequency for input 
to the converter, the 110-volt lighting circuit at the emergency generator 
was supplied from the emergency board. This total lighting load is less 
than five KW so that the effect on the alternator is negligible. The al- 
ternator load is assumed to from the bow thruster motor only. 

TERMINAL VOLTAGE The terminal voltage ae recorded using a Clevite 
Corporation Volts Converter. The recorder delivers a D.C. voltage propor- 
tional to the peak A.C. voltage. The converter was driven by the 220 volt 
output of the load sensing network potential transformer. 

ARMATURE CURRENT The armature current E recorded from the voltage 
across the input resistors to the load sensing network. The current to 
these resistors is supplied by the current transformers connected across 
the main line. Since this signal was at the 60 Hz line frequency the 
maximum amplitude of the recorded signal was limited to approximately 
twenty percent of the chart deflection in order to have satisfactory fre- 
quency response. 

BOW THRUSTER SPEED The induction motor drives an A.C. tachometer 
generator. The tachometer output was recorded directly on the Brush Re- 
corder. Speed can be determined either from the instantaneous frequency 


or voltage amplitude. For 1200 motor RPM the tachometer frequency is 





OG = 
40 Hz. The tachometer signal was recorded full scale on the chart since 


the recorder frequency response is flat to 40 Hz. 
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4.1 GAS TURBINE MODEL The derivation and final structure of the 
gas turbine model depends on the type of study involved. In one type of 
study, knowledge of internal behavior in response to MR. in environ- 
ment and load is required. In other EE knowledge of the time 
history of the input/output variables is required. Obviously for the 
latter case the model can be simpler and will most likely be in the form 
of a transfer function. 
In reference (10) a model is derived from basic principles for a 
Single shaft gas turbine where it is shown that the gas turbine can be 
represented by a multiple input-multiple output transfer function of the 


following form. 








5 
1 Gas T 
11 
Turbine 
Pu 
Engine 
N ら = 
m gb 


The matrix representing the input and output vectors is a relatively com- 
plicated non-linear function of the internal variables and engine charac- 
teristics. For this case the model can be simplified considerably and 
will be used in a linearized form. 

The most valuable information provided by EE e is an indi- 
cation of the functional relationships which exist between the input and 
output variables. Since the mass flow rate (m) is not being considered, 
only the functional relationship for torque is needed. From the transfer 
function it is seen that 


Ger 0 T1; Py, N, Mfp) 
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For this study, the ambient conditions P, and T, can be assumed constant. 
For a constant mass flow rate, fixed geometry and constant ambient condi- 
ES hé exit pressure Pr is also relatively constant. Since the mass 
flow rate is primarily a function of the speed, which is relatively con- 
Stant, the only input variable which changes significantly is the fuel flow 
rate me. As a result the functional relationship can be reduced to: 
て = $ (mrp) 
This functional relationship can also be derived by considering a 


control surface around the engine. 


Air Flow | Air Flow 


Engine 

Fuel Flow Shaft Power 
For this case where the engine can be considered small with respect to its 
ability to store energy as a result of either mass or thermal capacity, 
the changes in stored energy can be EN From observations on the 
prototype, it is known that the temperature of the exhaust gas is relative- 
ly constant throughout the operating range. With a constant air flow rate 
and exit velocity the enthalpy of the entering and leaving gases is rela- 
tively constant. As a result, the primary energy quantities entering and 
leaving the control surface are in the fuel flow and shaft work. With all 
of the assumptions above it is seen that output power is proportional to 
input fuel. Since speed is relatively constant it follows that 

T= (m) = Kie 
This relationship is verified by the engine performance curves given 


in Figure 4. Since the output KW is proportional to torque at constant 
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speed, the torque is proportional to fuel flow rate for constant Ti. 

In a simulation, it is necessary to consider the effects of the 
engine dead time or transport delay. This delay is defined as the time 
interval between a change in fuel flow rate and the change in torque out- 
put. From the control surface above, it is seen that this delay depends 
on the ability of the engine to store energy and the speed at which changes 
in internal temperature resulting from fuel rate changes can propogate 
through the engine. As stated above, the storage effects can be neglected. 
An estimate of the delay due to propagation time can be made by consider- 
ing the transit time for sound waves through the engine. This is the speed 
at which the temperature changes resulting from fuel rate changes propogate. 
For the Saturn engine, the transit or dead time is on the order of 0.005 


seconds. The torque relationship is then modified as follows: 


K! m 
T muse ~ ST “fb 
T = Ke Mp 7 + Ts 


The delay is approximated by a first order lag for the purposes of analog 
computation. 

All that is necessary to develop the turbine model is to determine the 
torque coefficient from the performance curve. This coefficient is a func- 
tion of ambient temperature with K, decreasing with increasing ambient tem- 
perature. It will be shown below that when the system block diagram is 
Mee timensicnalized for the analog computation, the torque constant can be 
lumped with the overall loop gain. 

4.2 GOVERNOR SECTION The governor schematic is shown in Figure 5. 
The engine block diagram can be formulated directly from the schematic us- 


ing the governor transfer functions supplied by Woodward Governor Company. 
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Figure 12. 


From the diagram it 1s seen that if the delay is neglected 

D oz SK NEO K 'Kr 9 
where it is assumed that the fuel flow is proportional to the angle of the 
fuel control valve. Although in most cases this relationship is approxi- 
mately linear, slight non-linearities will not affect the overall system 
performance. Although the fuel control valve is actually a dynamic opera- 
tor, in this case the valve will be considered to be a static operator. 
For low frequency changes in fuel angle this is a good approximation. 
With K, = K,'Ke ; res K,O Taking て = rsrr when O = (Omax - Omin)> 
where (Omax - Omin? is the maximum travel of the fuel control actuator, 
the torque constant K, is 


_ TREF 
Omax ^ Omin 
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In reference (10) it is shown that if the loop gain K, is not exces- 


sive the gain and phase margins at the crossover point of the open loop 


transfer function (approximately 6 rad/sec) will be adequate and the action 


1 


— can be neglected. 
SP i ye و‎ 


of the stabilizer section of the control box ( 
This simplification neglects the response of the governor system to fre- 
quencies above approximately 50 rad/sec. With the stabilizer section ne- 


glected and using the nominal values for the governor settings, the system 


block diagram is as shown in the section on Analog Scaling page 58. 


4.3 ALTERNATOR MODEL -The study of alternators typically involves 
deriving an equivalent circuit and phasor diagram of the electrical varia- 
bles. However, the phasor diagram does not explicitly show the mechanical- 
electrical energy conversion process which is a major part of the prototype 
system's power interactions. The model used for the simulation must in- 
clude the torque-current and speed-voltage relationships involved in the 
energy conversion process. Such a model can be derived from a study of 
the energy interactions of idealized magnetic circuits. The model proposed 
by Paynter, reference (10),was formulated this way and as a result shows 
the structure of alternator as a gyrator-type converter as well as the 
behavior of the electrical variables. 

The major assumptions required to simplify the analytical derivation 
of the alternator model are that the magnetic circuits are not saturated 
and that all fluxes are linear functions of abe currents. The lineariza- 
tion reduces the problem to a study of linear-coupled circuits and allows 
superposition of the magnetic fields. It is also assumed that the air 


gap is uniform. However, in a study of terminal voltage, current, torque, 
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and speed relations under load conditions, salient pole effects are rel- 
atively minor and can be neglected. One of the main effects of saliency 
is to produce a reluctance torque, and, except at no load, this torque 
is negligible compared to the total torque. The resulting model can then 
be used for machines which are operated under quasi-linear conditions and 
can be used for modeling the prototype alternator. 
The basic variables and relations required for the derivation are as 

follows: 

à = flux linkage = Nọ 

N = effective turns ratio of winding 


equivalent flux linking all terms = PF 


تد 
III‏ 


P = permeance 


n 


magnetomotive force = Nj 


Wm“ 


ul 


co-energy in magnetic field - [ «(Das = Pana; 

The inherent behavior of polyphase machines is the same regardless of 
the number of phases. As a result the torque is determined by the current 
flowing in each phase and the number of phases. In the balanced po1yphase 
machine each phase torque is a pulsating torque at twice the line fre- 
quency with a D.C. level which is exactly the total torque divided by the 
number of phases. When all phase torques are summed the pulsating torques 
cancel so that the total torque is constant in the steady state. 

A reasonably accurate model of the alternator (and motor) can be made 
using only a single phase if the torque pulses can be neglected.  Alternators 
with large inertias behave dynamically as low pass filters as seen from the 


simplified bond graph and system equations for a rotating mass. 
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Where > Where 0-7 


Is + R £ 
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For the prototype system I is on the order of 2000 lb-ft“ and 2w 


the frequency of the torque pulsations is 754 radians per second. Since 
R is relatively small; 


where w is the fluctuation in speed 
Ex | sen 
> x 754 caused by the pulsating torque. 
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The inertia of the induction motor and load is the same order of mag- 
nitude so that a single phase model can be used for both with torque pulsa- 
tions being filtered by the relatively large inertias. 

The alternator model used below is taken from one phase of a two- 
phase alternator model. Although the model is derived for a two-phase 
machine, one phase can be used to represent one phase of a three-phase 
machine since the torque pulsations are both double-frequency pulsations 
and are filtered by the inertias. All that is required to use this model 
for a three-phase machine is to properly scale the D.C. level of the 
torque. 

The two-phase representation shown below is very similar to the Park's 
dqo transformation used in studying three-phase machines. In this case the 
variables are referred to the two-phase stator whereas in the dqo trans- 


formation the variables are viewed with respect to the rotor. 
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O 
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Total IU tanking coil a 
Ma = NaP[Ngi, + Nrigc} (1) 
Meal flux linking coil b 
Mp = NaP[Nait + NEirS] (2) 
EESTI flux linking coil f 
Ae 7 NgP[Nai4C + NaipS + Neie] (3) 
Where Na = Np M = NaNFP 
C = ) 9 N= Na? 
S =sin 8 


and the permeance is independent of O in accordance with the assumption 
of a smooth air gap. 

Now v“ = f[2adiz + Apdip + Agdif] 

Wnt = f{PNiadi, + CMigdig + PNip,di, + SMigdip + CMiadif + SMipdif 


I le Nas ebe + Mig[Ci, + Sip] (4) 
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For any electro-mechanical converter 


Wom = T so that s ۷)۷ وپ سس‎ 


The exact relationship between T, w, e and i is found from a con- 
sideration of power balance. 

Assuming a lossless system: 

Pmech = Pac (developed at air gap) 


P = Pac - Pstoreq (in magnetic field) 


elect output 
The alternator is in effect separated into two processes, energy conver- 
sion and energy storage. Equating power at the air gap yields the fol- 
lowing identity: 

TW = Pac 

Where for a two-phase machine 

Tw 7 eglaà t epiy 


Migo [-Si, + Cib] = esia + epi, 


or eg = [-SMi¢] w= Ra w 
ep = [CMig] o7 Ryo 
mre Ko > - SMif and Ry - CMig 


ım matrix notation 


ea a a 
ep = O O Ry "no 
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If < is separated into its two components the gyrator-like action of the 
energy conversion is evident. 
Ta = Ral e, = Raw 


a a 


Th = Rob eh = Rp w 
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The magnetic energy can be modeled as an I field. The three currents 
AE) and ir are the sources of the stored magnetic energy. Since 
ig << iz and 
if << e 
most of the energy is stored as a result of the armature currents so that 
changes in armature currents will account for most of the change in stored 


energy. This implies a model of the form 


E Ye 
la | la 
which is in fact the equivalent circuit of an alternator with negligible 


armature resistance. 


VFO = variable 
frequency 
oscillator 
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where E represents the prime mover. The causality of E is determined by‏ 
the structure of the complete system. The Analog simulation of this‏ 
model will represent the exact alternator behavior (within the limits of‏ 
original assumptions). However, the model requires six multipliers. The‏ 
need to minimize Analog hardware required for simulation makes it neces-‏ 
sary to use the single-phase model.‏ 

For this case 
pS Ral 


a a 


-MIr sin (wt) ia 


wt 


E s I. sin(wt -$) where O 


¢ = phase angle between i, and ej 


ei 
l! 


ーー MIgl, sin(wt) sin (ut -$) 


MI£I, sin (wt) sin (wt - $) 


“MIcIa [cos$ - cos(2 wt - $)] 


۳ئ٣٣٢‎ t + 1 = MI -I, COSY, 


T4 alone gives a D.C. value equal to % of the actual steady torque 
plus a double frequency pulsating value. As shown above the altemator 
rotating mass has a large amount of inertia and acts as a filter for the 
high frequency torque pulsations. As a result, the one phase approxima- 
tion can be used without any appreciable loss in accuracy. 

The total torque input to the prime mover can then be obtained merely 


by doubling the torque calculated above. 


For the simplified model the bond graph becomes: 
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In this model the inertia of the rotating mass has been included and is 
shown connected to the -l junction of the alternator. Now the causality 
required for the prime mover can be determined. 

The I element imposes the velocity at the 1 junction so that the fol- 


lowing causality must be imposed. 


eq 


E ww 

This requires that the prime mover be modeled as some type of torque 
source. This is as expected because if a velocity source were assumed, 
step changes in the input would impose a step change in the velocity of 

the rotating mass and to do so would require infinite power. As a result 


the model of the gas turbine must be of a form such that velocity is the 
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input signal and torque is the output signal. 
From the simplified model bond graph the state variable representa- 


tion can be formulated from the following equations. 


dus _ 1 کہ‎ 
EC uu s TD e = Ry 
eq " 

Ra 

dr: 1 ( 

po — (e-v) 

dt JE: R = MIf cos wt 
v = Zi 


This model involves five parameters. 


Teg = inertia of rotating masses 
I, s synchronous reactance of alternator 
M = this parameter is defined above as Dale owever, it is 


closely approximated by the slope of the air gap voltage 
line of the alternator at synchronous speed. 


(cos ut o 


EC w = constant 
| M 


lr 
661 01207:7607 for the initial simulation will be assumed 
constant. 


の connected load. 


The Analog setup sheet is shown on page 57 where the load Z is re- 


placed by the induction motor model. 
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4.4 INDUCTION MOTOR A commonly used induction motor equivalent 


circuit is the Steinmetz equivalent circuit shown below. 


R っ /S 





The derivation of this equivalent circuit is given in standard Electrical 
Engineering texts. This model can be simplified considerably without 
Bros error by neglecting the magnetizing reactance X,. The circuit 
can then be simplified to the following. 


Kao A ex 
Kee eq 1 2 
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R2 
٦ ES and Ry へ Ro 
し 
Except at standstill Req % Ro/s. In the case of the wound rotor induc- 
tion motor large resistance is added to the rotor circuit for starting so 
that through the entire operating range the equivalent resistance is 
Ro! | 
Req ^v 7A where R, = Ro + Raaded- 

For the case of the prototype motor, the added resistance is approxi- 
mately thirty times as large as R, and is shunted out in four steps so 
EB can be neglected at all speeds. 

By applying a power balance across the air gap and neglecting losses 
an expression for the developed torque can be derived. 

T = qI2? Ro! 

m Ws S 
where (w,) is the synchronous mechanical speed and (q) is the number of 


phases. 
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4,5 PROPELLER The torque relationship for any non-cavitating 
propeller can be represented by an equation of the form 
DE pin: c_ (B) 
where the water density (p) and diameter (D) are constants. The torque 
coefficient C,(8) is a function of the velocity angle (R) where B is de- 


fined as ß = tan T (V/V, ) 


where Vo is the axial velocity of flow through 
the propeller and 

Vw is the tangential velocity of the propeller at the maximum diameter. 
The angle B is identical to the propeller advance ratio (J) used in propel- 
ler design by naval ers where J = Vo/mnD. With B constant the 
propeller can be modeled by an expression E the form 

て -= inr | 

From the picture of the propeller shown in Figure 7, it is seen that 
the propeller operates in a nozzle with an axial length approximately one- 
half diameter. With this nozzle the flow rate through the propeller is 
of the approximate form Vo «Q/D* 

D cum 

As a result the velocity angle ËB becomes 


B = tan (FZ) = const 


Tn D 
The torque expression given above then reduces to the form commonly 
used to represent propeller loads over a wide range of RPM where torque 
is taken as directly proportional to RPM squared. 
The prototype system was operated with the ship tied to the pier 


when the data was taken. Although there is some slight motion of the 


ship during the operation of the bow thruster it was assumed that the 
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propeller is stationary and that the axial velocity of water through the 
propeller is due to water motion only. 

Peo BEARING RESISTANCES All bearing resistances will be neglected 
in the simulation. In the case of the alternator where the speed devia- 
tions are slight the bearing drag is constant and will have no significant 
influence. For the induction motor the resistance will vary, but the 
resistance effects are assumed to be small compared to the load and driving 


torques. The resistances can therefore be neglected with no serious error. 
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The system model to be used for this 


simulation requires values for the following parameters: 


NE - engine inertia 

In - accessory drive inertia 

is - speed reducer inertia 

ite alternator inertia 

Xsa > Synchronous reactance of the alternator 

Xem ^ Synchronous reactance of the induction motor 
K rc on resistance of motor 

Pep added series rotor resistances 

T. - motor inertia 

Ewe engine torque/fuel flow ratio 


SES Ve nor parameters 


Where possible these parameters were obtained from information sup- 


plied by the equipment manufacturers. However, in some cases it was 


necessary to estimate parameters. 


I} The engine inertia is divided into three parts, the compressor, 


turbine and accessory drive package. Estimates for compressor and turbine 


mmem@tias were obtained from Solar. 


Compressor 
Turbine Stage (1) 
Turbine Stage (2) 


Turbine Stage (3) 


1.750 in-1b-sec< 
532 
‚564 


916 





3.762 in-lb-sec? at 22,300 RPM 
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I2 The accessory drive is approximated as follows. The idler gear 

and output gear are the largest rotating parts at the intermediate speed. 

These are approximated as steel discs one-half inch thick and eight inches 


in diameter. For a single disc, 


I: mr? 
TD 2Tp 
UL o wiere [or steel = 487 Jb/ft 
m = . T(6H) 487 __ - 0,0184 1b-sec</in 
u(2) 1728-384 
_ 1 2 : 2 
NE = 0.147 sun sec 


For two discs and considering the remaining parts as having the same inertia 
EN 0ngle disc, the total inertia of the accessory drive is approximately 
0.5 in-lb-sec? at 6010 RPM. 

E The speed reducer inertia is estimated by assuming that all 
inertias are negligible except that of the first-stage ring gear. Since 
all inertias will be referred to the output speed, both speed ratios and 
WR? values must be considered. The first-stage ring gear is approximated 
by a one-quarter inch steel disc eleven inches in diameter plus a ring two 
inches wide by one-quarter inch thick and eleven inches in diameter. The 
resulting inertia is approximately 1.0 in-lb-sec?. Although the inertias 
of the second-stage parts are of the same order of magnitude as the first 
stage ring gear, they become ہے‎ compared to the first stage 
parts when referred to the output speed. 

Zu The alternator rotor inertia was given on the manufacturer's 
drawings. 


WR2 = 951 lb-ft? or = 5 lb-ft-sec? 
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All inertias can now be referred to the output speed of 1200 RPM. 


Nturb 2 
ii = I —— 
lEQ 1 E es 
2 
3.76 22 6ء‎ 1| 2 
7 Ul | 1200 | = 108 lb-ft-sec 
E: N acc drive 
I, = I, eee te 
EQ N alt 


= 1 mal 


N 2 
- NE | 1s stage 
3EQ 3 Nalt 


1.0 (6010)? 
< -— x ` Et. 2 
12 m 2 lb-ft-sec 


È r 2 
Tung = Ty = 29.6 lb-ft-sec 


It is seen that the inertias of the engine and the alternator are the 
° ° . . 2 
dominant terms. The total inertia is taken -as 140 lb-ft-sec . 
Xsa The alternator synchronous reactance is given cn General 


Electric drawings as 1.19 ohms. For a line frequency of 60 Hz; 


L = pm 0.00316h 


Xem, Ra The motor reactance can be estimated from the equivalent 


circuit model used for the motor and the motor performance curves. 


Asm ANE |. S 
RL = — و‎ Ss = 
5 S S Ns 





From performance curves at the motor rating point: 


N 001 ۳ 


I = 440 amps s = 1.66% 


E = 440 volts 
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For the equivalent circuit used, the magnitudes of the phasor quantities 


are related by: 


E E 2 25 
ت0‎ 7 where Z2 = (R, EX ) 
Ke = Qi la Om 
E 
ے‎ = I^ l ohm and L = 0.00188h 


The rotor resistance is given in the General Electric drawings as 0.0117 ohms. 

In No figures were available in the manufacturer's drawings for the 
motor inertia. The inertia was estimated from a drawing of the rotor where 
the mass was approximated by a steel right circular cylinder fifteen inches 
long by twenty inches in diameter. The inertia estimate is 14.1 lb-ft-sec?. 
This is approximately fifty percent the inertia of the alternator and is a 
reasonable estimate. 

The propeller is approximated by a steel disc four feet in diameter 
and one inch thick. This gives a value of 30 lb-ft-sec?. When this iner- 
tia is referred to the motor speed through a speed ratio of approximately 
2:1 the equivalent inertia of the propeller becomes 7.5 1b-ft-sec?. 

E The ratio of steady state engine torque to fuel flow is obtained 
from the engine performance curves. For constant speed and constant ambient 
conditions the torque is linear with fuel flow throughout the operating 
range. In this study where fuel input is not observed it is not necessary 
to isolate this constant from the overall gain in the engine control loop. 
In sealing of the simulation it is shown that where torque and fuel actu- 


ator angle are considered in non-dimensional form this constant is reduced 


ES Unity. 
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St No information on the propeller was available from the manu- 
facturer. However, for the assumed propeller model which 15 
tp = کپ‎ 
a torque coefficient can be fit from the observed propeller performance. 
Using the maximum observed propeller speed for fitting the coefficient 
gives the following: Nnay 7 1125 RPM 


Neglecting frictional losses the motor torque at this speed is 1370 1b-ft 


so that 


T 
"P 
KD Soe 


0.00108 lb-ft/RPM? 

GOVERNOR SETTINGS There are five parameters connected with the 
engine governing system. These parameters can vary over a wide range 
depending on the governor settings. In order to determine the exact 
settings it would be necessary to measure the open loop frequency re- 
sponse of the control system. The ranges of the governor settings as 
well as nominal values were supplied by the Woodward Governor Company. 
The nominal values will be used in this study. The ۳7ع 8۶ف‎ gain and 
loop gain will be varied and the effect on system performance observed. 


The values for the other parameters will be taken as follows: 


Parameter Range Nominal Value 
01 23 to 184 28 
a2 7 6 
ag 0.08 to 160 4 
K vary in simulation 305000 
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Parameter Range Nominal Value 
Y. == 46 
O ーー 2210 
Feed Forward Gain vary in simulation 0.0045 
Ra_ The added rotor resistances are given in the manufacturer's 


drawing. The resistances are shunted out in four steps. 


Step 1 .240 ohms/phase 
Step 2 .088 ohms/phase 
Step 3 .072 ohms/phase 


Step 4 .063 ohms/phase 
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5.2 SIMULATION-ANALOG SET UP The models derived in previous sec- 
tions can be set up for Analog computation in a straightforward manner. 
The set up for the complete system is shown in Figure 9 on page 57. 

GAS TURBINE The Analog set up for the gas turbine and speed control 
is implemented directly from the block diagram of the engine on page 31. 

(1) The feedforward of load disturbance is a high pass filter of the 


form: 


KS 
)و0۰‎ 


E 
X 
This transfer function is realized by a circuit of the form below 


which is taken from reference (3). 


٦ 


where A 


DAT 





The operator is realized using ABO and A31 in the simulation. 
(2) The engine speed is obtained at A60 by integrating the difference 


of the driving and load torques. 
ES 
Au = I | Gis _ Es Gt = wet, ) zu (to) 
to 
Since w (to) can be taken as the base speed of the alternator the 
output of A60 when properly scaled gives the error in engine speed. This 
voltage is summed with the output of the feedforward network to give the 


error signal developed in the EGA control box. 


(3) The error signal is generated by A61. 


K4 s m A = 【 
Ei = .. Ta “Ky に where T A 
ri Ure f a ret 
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(4) Since the operator representing the stabilizer section of the 
control box has been neglected the transfer function between the error 


signal E, and the output of the fuel actuator.angle O is 


© _K, (Tjs + 1) (Tas + 1) 
Ex (TE 1) 


This transfer function can be realized by taking two operators cas- 


caded. 
T +1 
For ae جا ری‎ 
x S 
I TX + X 
m Y = T,X + JX 


L 
This operator is realized using A95, A39, and A01 





Tel) 
The operator for x^ —£ — is also from reference (3). 


(T,s + 1) 
A = Klo/Tz 
B = (T。-T。)/T2T。 
C = 1/T4 





The output of A90 represents the angle of the fuel control actuator. In 
the derivation of the turbine model it was shown that except for a trans- 


port delay the torque is proportional to fuel input. The fuel control 





eg 

angle is assumed proportional to fuel input. The engine dead time can 
be neglected so that except for a gain factor the output of A90 represents 
the gas turbine tcrque. 

With the correct setting on P60 the e c of A91 represents engine 
torque. 

The gas turbine model is now complete except for the load torque. The 
alternator torque is generated bv the multiplier A38. This torque is of 
the form 

t 7 t9 (1 * cos 2 ut) 
It is necessary to filter the double frequency component for two reasons. 

(1) The pulsations are a form of noise which cause slight variations 
in engine speed. When compared to the base engine speed of 1200 RPM, 
these fluctuations in speed are insignificant. However, with a fine speed 
control on the engine where the variations in speed due to the real load 
changes are of the order of 10 RPM, the noise from the load torque signal 
becomes significant. This point is most evident when scaling the integrator 
which generates the speed signal. If the integrator is scaled so that the 
output is Dk ۰.۱۹۷000 EEN the integrator gain which represents 
the reciprocal of inertia is small and corresponds to a large inertia. In 
the simulation it is necessary to generate N/N' where N' is of the order of 
the maximum deviation in engine speed due to changes in load torque. As a 
result, the integrator gain is increased to the point where the double 
frequency torque pulsations are passed along with the D.C. torque level. 

(2) It is also necessary to filter the load torque because it is in- 


put to the load sensing network which has a derivative action. 
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A low pass filter was chosen over a band stop filter for removing 
the torque pulsations because the low pass filter requires less Analog 
hardware. An attenuation of 40 db at the noise frequency will provide 
adequate filtering. At this level the noise is barely noticeable on 
the oscilloscope. A first order filter tuned for an attenuation of 40 
db at the noise frequency of 120 Hz will have its break point at 1.2 Hz 
anc is considered to be too low for the expected load transients. As a 
result a second order filter was used with the break point at 12 Hz, and 
40 db attenuation at 120 Hz. The transfer function of the second order 


filter is of the following form: 
1 
1 2 TI | ie 2 
On un) 


With the break point at 12 Hz and with Z set at 0.55 in order to minimize 


G(jw) = 


overshoot near the break frequency, the transfer function becomes: 


1 
G(s) = 


> | = 


> where s - jw 
5 Š رد کھ ا‎ 
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Performing the operations indicated gives 
Y = 5650X - 82.5 Y -5650Y. 


Scaling this equation to a form suitable for Analog computation gives 


Y = SDA 3205 ے‎ Oc 
0 = : E : 100 = +7 


This filter is shown on the Analog set up sheet using A70, A74 and Al00. 
ALTERNATOR, INDUCTION MOTOR, PROPELLER The equations representing 


these components were derived previously. 





Ze 


er 
= 
| 


ge), = 5 Cee ts) 


C) R = oos we 


MS Ta = R 


e v) 


(5) e = Ru 
i^ i |Fouo iR N 
(6) v = — = —— = 6ت‎ 
S WO —W l. 


E 
7) pe ` ES Een 


(8) t = KN? 


P 
dN 30 
(9) Ee > 0 Ti OD ー 15! 


Equation (1) has been described above with A60 generating the 
alternator speed. | 

Equation (2) determines the gyrator modulus of the alternator and 
is a function of the field current and frequency. In the prototype the 
speed variations were insignificant compared to the base speed so that 
the frequency can be considered constant. Although the field current 
did vary slightly a constant field current will be used. The gyrator 
modulus is then of the form 

R = K cos (wt) 
The 60 Hz sinusoidal waveform was generated using the Wavetek signal 
generator. 

Equation (4) is solved using AOO where e is the air gap voltage of 


Aia ltermator and v is the back emf of the induction motor. 





EN 

Equation (5) is the air gap voltage of the alternator. With the 
assumption of constant frequency and field C M e is a sinusoidal 
voltage of constant frequency and amplitude. 

Equation (6) generates the back emf of the induction motor and is 
realized at multiplier AO3. 

Equation (7) generates the torque output of the induction motor. 
This torque also contains a double frequency component. However, since 
the induction motor speed is scaled to the 1200 RPM synchronous speed 
it is not necessary to filter this torque. 

Equation (8) is realized at multiplier A63 where the propeller 
torque coefficient is set on P67. 

The switching of the induction motor rotor resistances is set up 
in the network which is input to the junction inverter AO4. The values 
of the switched resistors are set on the potentiometers P95, P00, P32, 
and P65 with the residual resistance on P03. The switched resistors are 
input to the amplifier junction of AOà through the D/A switches shown. 
The switching of the D/A switches is controlled by the logic outputs of 
the comparators. 

In the prototype system the switching is controlled by time delay 
relays. This switching action is simulated in the following way. 05ھ‎ 
generates a ramp timing signal. This signal is input to comparators C04, 
C34, C65, and C94. These comparators have a negative bias corresponding 
to the time at which switching is to occur. The inverted logic output 
of each comparator is input to the switch control of the respective D/A 


switch. Initially the output of A05 is zero. All comparators are in the 





EDS 
low position. As a result, the inverted output to the D/A switch is high 
and all switches are conducting. The input to AO4 is 

IR + OR] + Ro + Ra + Ry). 
As the output of A05 increases to a level ne the bias of each comparator, 
switching occurs when the comparator logic output goes high and the in- 
verted output goes low so that the respective D/A switches are turned off. 
After all switching occurs, only the resistance due to PO3 remains in the 


ei euit. 
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5.3 SCALING ANALOG SIMULATION 
GAS TURBINE When nominal values are substituted for the govemor 
settings the block diagram for the prime mover and governor is as shown 
below. 


M ۲۲076 13. 


0.00405s 


(s/10+1) 





1 1035(s/10 + 1)(s/4 + 1) 
» Tp 5 ü ^ s(s/28 * 1) 





En 


The machine time constant EE is defined from the non-dimensional equa- 


tion for the alternator speed. 


ET Iw 
le "a = wo Es I where T - = 
To TO dt wo d 0 


The base speed Du is 125.6 rad/sec and the inertia I is the combined iner- 
Eta Of the engine, speed reducer and alternator referred to w,. The selec- 
tion of the base torque is somewhat arbitrary. Based on the engine power 


rating of 1200 HP with output at 1200 RPM the base torque is: 


000 0 
To see ٠× 5250 1b-ft 


The maximum expected steady load torque is approximately 1500 lb-ft. The 





e 
maximum engine torque is limited by the acceleration limiting response 
of the fuel control valve. For a step load change of 1500 lb-ft the 
maximum engine torque developed is assumed to be 3000 1b-ft. With this 
base torque the engine time constant = 9.86 see. The EE gain 
on the integrator which generates Au is 0.170. The expected maximum 
deviation in speed is od be 7.5 rad/sec. With this base speed the 


gain on the speed integrator is: 


Aw 1 Ta T 





. = = a 
85 Tm's o 
X 
where T,' - . < یی‎ me 


The gain on the speed integrator must be increased by a factor of 16.7 
for this new value of of Tm'. This is accomplished by modifying the sys- 


tem transfer function as shown below. 








0.0450s 
(s/10 * 1) 







CNET) 
s(s/28 + 1) 
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The operator Tt. is generated by A60. The scaled equation is* 
m 


Ao | _ 10(2840) |Ta - Ta 
05 S 3000 
with the pot settings on P38 and P97. 


The scaled feedforward transfer function with T = 10 and K = 0.1 is 


generated at A40 and A31. The gain of the feedforward is then set on 007. 


(s/4 + 1) 
5 


The operator is generated at AOL. For T, = 4 the setting 


1 
۲۰۲١۹ is 2500. 
×۴ ۶۹/10 + 1( 
The operator [EE generated at A90. In order to reduce 
the number of pots required K is taken as 100T3/T,. With T, = 1/10 and 


Ta = 1/28 the pot settings and gains are: 


62 


29002800) 


7 


100(1800) 

The inputs to A65 and A91 from AOl are then at a gain of 100. Since the 
gain of 100 is difficult to obtain at the کس‎ to summer A91 the output 
of A91 is scaled down by a factor of 10. The inputs and output of A91 


are then as shown on the Analog diagram, Figure 9. The output of A91 cor- 


1 
responds to 1o ES 
12 
Hie Loop gain is introduced at A90. With K = ج‎ and 
max min 


Onax - Onin 


taken as 30° the required gain at P60 is determined as shown 
* The following convention is used to represent all scaled equations: 
normalized variables in machine units are in brackets, potentiometer co- 


efficients are in parentheses and amplifier gains before the parenthesis. 
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below. 
92.5 KO _ 92.50 
Te Eo Omax 7 Omin 


The output of A91 is 10 2 € 


2 
mete (gain) Oo = 92.5 0 
To Omax 7min 
(gain) = 2 


— - ی٤‎ 
T3( Omax ü Onin? 


The required setting on P60 is then 1(8650). 
LOAD COMPONENTS The induction motor speed is generated at A35. 


The unscaled equation is: 


The maximum motor speed is 1200 RPM which is the synchronous speed. The 
maximum steady motor torque is approximately 1500 lb-ft. Since the simu- 
lated motor torque has the double frequency component the maximum torque 
is twice the instantaneous D.C. level. As a result the maximum torque 
is taken as 4000 lb-ft. With the inertia (I) taken as 21.6 lb-ft-sec* 


the scaled equation becomes: 


DNO (3875). | im = Tp 
1200 5 s 4000 Set on oo 


The propeller torque is generated at A63. The unscaled equation is: 
= 2 
تم‎ KpN 
With E - 1.08 x 10^ and the torque and speed base values as above the 


scaled equation becomes: 


ES :(3880) N< 
4000 ۰007 756 Set on P67 


سے 
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Ihe voltages and current in the simulation must be scaled based on 
peak values rather than the RMS values. For maximum RMS values of 440 


volts and I/phase of 1000 amps the base values are taken as: 


Ese 1414 amps 

Vbase = 750 volts 

Multiplier A33 generates the unscaled product of e and i. This is 
scaled to the induction motor torque as follows. For the induction motor 
for three phases: 


met) = 3 V SG) sie sum 71+ 


Max 
or 


3 V I (1 + cos?wt) 
See ele le eee eee 
Belt) dm 


where the phase anzle can be neglected in scaling the equation and (.7376) 
femeconversion from newton-meters to lb-ft. w, is the synchronous me- 
chanical speed which for this case is 125.6 rad/sec. With the base vari- 


ables assumed above the scaled equation becomes: 


ES 2 30(2310)|], V. | L1. 
4000] ` 7501 | 151% Set ON LG 


The output of A00 must be scaled to Imax 7 1414. The unscaled equa- 


Lion for A00 is: 


à ] (e -V ) 
Jr جو‎ ecce سا‎ 
S 


L 
Using the base values for current and voltage and the equivalent inductance 


(L) for a single phase the scaled equation is: 


i} _ 10° (8460) | e -v 
1414 3 750 


The motor voltage is generated at A03. The unscaled equation is: 





RaINs _ RI 
= Ns-N ^ 1-N/N。 








ge 
The denominato» is generated at A06 and is already in a non-dimensional 


form. Scaling the equation for current and voltage gives: 


E 1.89 R4 i 
750 S y 1414 


The potentiometers in the switching network are then set to 1.89 R, and 
Ra' where Ra is the value for the rotor resistance and Ra' the value of 
the added rotor resistance. 

The alternator torque is generated at multiplier A38. As shown 
above this output can be scaled to a base torque of 4000 lb-ft using the 
same setting as on P36 SEDE 10(2340). However, the alternator torque 
is filtered before being input to the gas turbine section so that it is 
not necessary to scale this signal for the EE frequency pulsation. With 
the turbine base torque taken as 3000 lb-ft the gain on the output of A38 
MIS 100). Since the second order filter is set to a gain of (1) this 
gain is set on ٠ 

The scaled variables for all amplifiers and settings for all poten- 


tiometers are given on the assignment sheets given in Tables 1 and 2. 





Amplifier 


00 
01 
02 
03 
oy 


05 


06 
07 
09 
10 
14 
30 
31 
32 
33 
34 
35 
36 
38 
39 


40 


AMPLIFTER ASSIGNMENTS 


Scaled Value 


7/10 


e 7750 
v/750 


7 +0 


(N。-N)/N 
2773090 
v/750 
(v-e)/750 


-(N。-N)/N 


-i/1414 
vi/750 x 1414 
-N7/1.44 x 108 
N/1200 


(€ e 075 


P 


Comment 


Governor Control Box 
Alternator Modulus 
I.M. Back emf 

For I.M. Back emf 


Time Signal for Resistor 
Network 


5۰٠۰16 


Invert I.M. Back emf 


Invert A06 


LPF for Viewing Torques 


FFWD of Governor 


I.M. Torque 


I.M. Speed 


Load Torque on G.T. 


Governor Control Box 


Governor FFWD Net 





Amplifier 


44 


60 


61 


63 


65 


70 


74 


20 


BEI 


25 


100 


Scaled Value 


Aw/7.5 


N2/1.44 x 108 


+。/3000 


ES 


Comment 
Invert A02 
Change in Alternator Speed 
Governor Error Signal 
For Propeller Torque 
Governor Control Box 
Butterworth Filter 
Butterworth Filter 
Es E 0)6 
Fuel Control Angle 
Governor Control Box 


Butterworth Filter 





EOS 


Table 2. 


POTENTIOMETER ASSIGNMENTS 


POT PARAMETER SETTING W/GAIN 
QOH Aw/wg Gain 1000 (10) 

007 FFWD gain 2000 (vary) 

00 Ra (switched) 6000 

02 RIE. modulus cet fOr Ee 
03 Ra 0597 

04 Switching time 1002 

05 Timing signal 0500 

07 Scale Ta +3 ۳ 

08 Cont. box setting 2287 

10 Motor inductance 8277 (1000) 

32 Ra (switched) 2000 

34 Switching time 3291 

35 Inertia of I.M. 2000 (10) 

36 Scale motor torque و‎ 210) 

ey Gove Cont. box 1800 (10) 

38 Alt. inertia 2841 (10) 

60 Control loop ٤0 2177 (10) (vary) 
61 Butterworth filter 5649 (100) 

62 Govt. Cont. box 286120100) 

64 Switching time 4895 

65 Ra (switched) 2005 

E Scale propeller torque 3880 

70 Butterworth filter 8243 

90 Switching time 6819 

95 Ra (switched) 2001 

97 Alt. inertia 2841 (10) 

100 Butterworth filter 5648 (100) 
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CG. RESULTS 


The results of monitoring the prototype performance and the results 
of the simulation are shown on pages 69 through 75. A chart speed of 
twenty-five millimeters/second was used for all recordings. Strips A-L 
correspond to a time duration of approximately fourteen seconds while 
channels M-BB span approximately the first eight seconds of the simulation. 
The records are presented in the followins order. 

A. Observed Alternator Speed Deviation 

B. Simulated Alternator Speed Deviation 

C. Observed K.W. Load (through 0.068 sec. LPF) 


D. Simulated Alternator Torque (through 0.068 sec. LPF) 


E. Observed Line Current 
F. Simulated Line Current 
G. Observed Induction Motor Speed (A.C. tachometer) 


H. Simulated Induction Motor Speed 


I. Simulated Motor Torque (through 0.068 sec. LPF) 
J. Simulated Alternator Torque (through Butterworth filter) 
K. Simulated Engine Torque - No Dead Time, No FFWD 


L. Simulated Engine Torque - No Dead Time, With FFWD 


M. Alternator Speed with FFWD Gain (Q07) at .1000 
N. Alternator Speed with FFWD Gain at .2000 
O. Alternator Speed with FFWD Gain at .3000 


P. Alternator Speed with FFWD Gain at .5000 





mc 


DD. 


Simulated 


Simulated 


error 
Output of Governor Load Sensing Section 


Governor Error Signal 


Observed Alternator Terminal Voltage 


Simulated 


Simulated 


Simulated 


Simulated 


Simulated 


Simulated 


Simulated 


Simulated 


Alternator Speed with 0.1 Sec. Dead Time 
Engine Torque with 0.1 Sec. Dead Time 
Fuel Control Angle with 0.1 Sec. Dead Time 


Alternator Speed with 0.01 Sec. Dead Time 


Engine Torque with 0.01 Sec. Dead Time 
Alternator Speed - Loop Gain = 1.0 Per Unit 
Alternator Speed - Loop Gain = 1.6 Per Unit 


Alternator Speed - Loop Gain = 2.3 Per Unit 
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DISCUSSION OF RESULTS 


ALTERNATOR SPEED A comparison of strips A and B shows that there is 
excellent qualitative correlation between the simulation and the prototype 
performance. In the simulation the speed transient which occurs when the 
induction motor switches to speed five occurs approximately one second too 
late. This is due to an error in setting P90 which controls the switching 
time for motor speed five and in no way affects the validity of the model. 
The three criteria listed below give an indication of the quantitative cor- 
relation for the first speed transient. 

(1) Maximum deviation in speed--For the prototype the drop in fre- 
quency is eighteen chart divisions at a sensitivity of 20 mv/division. The 
frequency converter sensitivity is 2.5V/5Hz. The frequency deviation is 
0.72 Hz which corresponds to a 0.24 rev/sec drop in alternator speed. In 
the simulation the speed drop was 1.5 rad/sec or 0.239 rev/sec. 

(2) Ratio of speed drop to speed overshoot--Comparing the simulation 
and prototype gives the following approximate values: 

Simulation - 3.34 
Erototype = 3.60 

(3) Time to return to original speed--For both simulation and proto- 
type the speed first returned to synchronous speed in approximately 0.70 
seconds. 

The Be elation noted above is admittedly a retrofit to the prototype 
data. However, it is clearly shown that the model used closely predicts 
the prototype performance. In the simulation the engine parameters set 


on the potentiometers were the values determined in the previous section. 
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The only settings which varied were the governor loop gain and the gain 
on the governor feedforward signal. These En settings are varied in the 
simulations to show their effects on engine performance. 

The simulated output of the feedforward network is shown on Seri pace 

As expected this signal shows the derivative operation in the FFWD net- 
work. Strips M through P show the effect of increasing the FFWD gain. 
Bae strip P the speed shows an initial increase when the load is applied. 
This type of response is not possible for the system under study so that 
an upper limit must be imposed on Q07 which corresponds to the FFWD gain. 
However, the trend in engine performance as Q07 is increased is evident. 
iMefettect of the FFWD on engine torque is pet in strips K and L. 

The effect of loop gain on performance is shown in strips Z, AA and 
BB. As gain is increased the speed deviations become smaller. In strips 
AA and BB it is evident that the FFWD gain is set too high since the speed 
Shows a slight initial increase as load is s However, the influence 
of loop gain is still clear. 

In deriving the turbine model, the engine dead time was assumed to be 
negligible or at least less than 0.01 seconds. The influence of dead time 
on response is shown in strips U through Y. The dead time is simulated us- 
ing a first order lag approximation. In strip U the speed response for 0.1 
seconds dead time is shown. The response is noticeably worse than that 
shown on strip N which is the response for the same potentiometer settings 
for all components and with no dead time. The reason for the deterioration 
in response is evident from comparing strips V and W. The dead time causes 
the torque to lag the fuel actuator angle as shown. Strips X and Y u 


Spond to a dead time of 0.01 second. Comparison of the speed response 
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shown on strips X and N shows that there is no noticeable difference in 
response. As a result, any dead time less than 0.01 second can be neglected. 

ALTERNATOR LOAD Strip C is the alternator KW load as measured by the 
governor load sensing section and as viewed through the 0.068 sec. LPF. 
Since the maximum alternator speed variation is on the order of one percent 
of base speed the load is roughly proportional to the alternator torque. 
The load signal was not calibrated to a known level of load so that this 
Measurement may be used only for a qualitative comparison. Two ratios are 
used below to give an indication of the degree of correlation. The over- 
shoot ratio compares the maximum torque for each transient to the steady 
torque level at the corresponding induction motor speed. The steady state 
ratio compares the steady torque at each induction motor speed to the steady 


torque at the maximum motor speed. 





Table 3. 
OVERSHOOT RATIO STEADY RATIO 
IM SPEED SIM PROTO SIM PROTO 
1 2.80 2.00 213 .55 
2 1.29 1 5 . 70 
3 1.16 30+0 12 .78 
Y: E 1.19 . SH 
M 230٦ 1.30 00 17.00 


The differences between simulated and observed values is due primarily to the 


difference between estimated and actual rotor resistance values. 


A rough check can be made on the actual magnitude of the simulated 
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alternator torque. At speed five the induction motor torque is approximate- 
ly 1500 lb-ft at 1150 RPM. Since the alternator speed is approximately 1200 
RPM and the power levels of the alternator and induction motor are approxi- 
mately equal, the alternator torque will be approximately 1400 lb-ft. 6 
scaled value for the simulated alternator torque is 1350 lb-ft. 

INDUCTION MOTOR The m of the induction motor performance 
and the simulation is given in Table H below. The values for the motor 
torque and current were taken from the motor performance curve, Figure 11. 


Also shown on this curve are the simulated values for torque and current. 


Table 4. 
Speed RPM RPM Torque Rs Torque LRMS 
sition Obs. Sim. (Performance Curve) Simulated 
J 834 840 800 = 700 200 
2 266 960 1000 S10 dus 250 
3 1082 1015 1100 340 1030 275 
L 1104 1080 1270 380 90 300 
o TEBO 1150 1450 425 1360 350 


As with the engine simulation the results can vary widely depending on the 
values used for the motor and propeller parameters. The original estimates 
for the motor inertia and the rotor resistances were changed slightly in 

order to provide a better fit between the model and the observed performance. 
The differences between the original parameter estimates and the values used 


in the simulation are summarized in Table 5. 
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Table 5. 
Estimated Setting From 
Parameter Value Potentiometer Scaled Ean. Setting Used 

a .71 ohms P10 8460(10?) 8277(10°) 
p DER P03 0220 0597 
Ra(1) . 240 POO 4450 6000 
Ra(2) .088 P32 1660 2000 
Ra(3) 2072 P65 1360 2000 

ESI .063 P95 1190 2000 
ImtI, 21.6 P35 | 7 0 2000(10) 
K 1.08 x 107 P67 3880 3880 


All parameters agree with the estimated values except for the rotor 
resistances and the motor inertia. The difference in inertia is not large 
and is within the accuracy of the original estimate. The potentiometer 
٣۰٠٠٦٠٠ corresponding to values of R4 were considerably higher than the 
original estimates. These higher settings result from the fact that a 
single-phase model is being used to represent a three-phase machine. 

For the three phase delta connected rotor the resistance is given as 
0.0117 ohms. This resistance was probably measured across the slip rings 


and would correspond to the following diagram. 





As a result the single phase equivalent resistance (Rp) ۹ 48 


requires that the settings for the potentiometers representing the rotor 
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resistances in the table above be increased by a factor of 1.5. The values 
given in the manufacturer's drawings for the switched resistors are -—À 
to be nominal values. For this simulation the same value (.70 ohms) was 
used for resistors two, three and four with resistor one assumed to be three 
times as large. These values are very close to the specified nominal values 
and give excellent qualitative correlation for the transient currents. 

With R, for the rotor given as 0.0117 ohms and corrected to a single 
phase equivalent the required setting on P03 is 0032 whereas a setting of 
0597 was used in the simulation. A value of 0032 was used but the current 
surge when switching to speed five was larger than observed on the prototype. 
As would be expected, the current transients at the other switching points 
were not noticeably affected by this change in the residual rotor resistance. 
The difference between the required rotor resistance potentiometer setting 
and the value actually used is due to the fact that at this point the divider 
AO3 which is generating the motor back emf is poorly scaled for Analog compu- 
tation and M acy is limited. 

When P65 is switched out of the motor resistance network the operation 


AOS 3s as follows: 


At the switching point (N, - N) is approximately 0.12 per unit and approaches 
0.05 P.U. in the limit. With (I) properly scaled to avoid saturation at all 
other transients the magnitude of (I) at switching is approximately 0.3 P.U. 
With R, set at 0300 the component A03 is dividing a small signal by a small 

Penal. As a result accuracy at this point is limited. A03 can be rescaled 


to give higher accuracy when N approaches N, but this will cause saturation 
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when the motor is first started. The divider must be constrained to 
operate at all times with the numerator smaller in magnitude than the 
denominator. 

Even with the reduced accuracy. when the motor speed approaches 
synchronous speed, there is still good correlation between the simulation 
and the prototype EE throughout the entire operating range. 
Strips E through H show the correlation of speed and current. The data 
tabulated above is plotted on the motor performance curve, Figure 11, and 
shows the quantitative correlation for torque and current. The simulated 
current shows the same trend as the performance curve data but is lower at 
all points. This is due in part to the error in measuring the magnitudes 


of the current from the strip recording. 
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7. CONCLUSIONS AND RECOMMENDATIONS 


The correlation between the simulation and prototype for all observed 
variables proves the validity of the models used in the simulation. The 
natural tendency is to expect that highly complex models will be required 
for each component. However, reasonable accuracy in quantitative predic- 
tions and highly valuable qualitative information about system performance 
have been obtained using simple "off the shelf" models. The simulation 
can be helpful in design or as a training aid. 

As indicated earlier, the simulation gives the designer an indication 
of trends in system performance as the various parameters are adjusted. For 
a given motor the parameters A cde fe o ii E sv tened resistors 
and switching times are the two most important parameters which can also be 
easily adjusted. The simulation shows how the performance can be varied 
over a wide range by adjusting these parameters. Although the designer has 
little control over the inertia of the motor and propeller, the inertia ef- 
fects are clearly shown. Finally, the effects of propeller pitch as char- 
acterized by the propeller torque coefficient can be investigated. 

The simulation provides an excellent study of the governing system 
and engine performance. The most important governor adjustments are con- 
trollable by potentiometer settings so that changes in the govemor coef- 
ficients can be investigated. The effects on engine torque and speed 
response are clearly shown. In addition e useful design informa- 
tion, the simulation can demonstrate to operating personnel the effects on 
system performance due to changes in the governor settings. 


The transfer function approach to modeling the prime mover yielded a 
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simulation which can easily be adapted to any other type of prime mover by 
substituting the proper transfer function Con a diode function generator 
if necessary) in the place of the gas turbine model. The model used in 
this simulation is readily adaptable for use with diesel engines with the 
most significant change being in the assumed engine dead time. 

The two most important dynamic effects which were not included in the 
simulation because of the limitation on computing hardware were the voltage 
regulation and fuel control valve effects. With additional hardward these 
effects can easily be included. However, as the results indicate, neglect 
of these effects did not seriously affect the accuracy of the simulation. 

Neglecting the stabilizer section of the governor control box and 
dynamics of the fuel control valve did not seriously affect the results 
for two reasons. The input disturbances to the gas turbine were below the 
frequency level where the stabilizer section influences performance. Al- 
though the load transient was close to a er change the second order filter 
on load inputs to the gas turbine limits change to approximately 12 Hz. Al- 
so, in scaling the simulation the developed engine torque was limited to 
3000 lb-ft. This provides a form of acceleration limiting. 

Neglect of the voltage regulation effects can influence the results 
in two ways. First, in the simulation the induction motor is supplied with 
a constant voltage whereas in the prototype the voltage varies as shown in 
Strip S. The error caused by this assumption is not serious. The second 
effect of neglecting the voltage regulation is on alternator torque. With 
field current assumed constant, the gyrator modulus on the alternator is 


constant so that the torque is directly proportional to the line current. 
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The field current does رر‎ sllghtly durinz the transients so that the 
prototype alternator torque is higher than shown on strip J. However, E 
correlation between the simulation and prototype performance shows that 
considering field current constant does not cause a serious error. 

The use of a single-phase model for the three-phase machine intro- 
duces no serious errors. The inertia of the induction motor is large 
enough so that the double frequency torque pulsations in the motor are not 
noticed. At the alternator end, the second order filter provided adequate 
filtering of the pulsations and adequate frequency response in order to pass 
the transient loads. Even with unlimited Analog hardware it would not be 
necessary to use a three-phase model as the single-phase model is more than. 


adequate. 
RECOMMENDATIONS 


The observed system performance indicates that the original system 
design is more than adequate in order to meet its performance specifica- 
tions. Therefore, use of this simulation in an attempt to improve the 
overall design is not practical. The only change to the original system 
which could be made at small cost is in the starting resistors. The speed 
range of the motor with the installed resistors is from approximately 65 
to 95 percent of synchronous speed. This range could be extended by chang- 
ing the resistors. The required resistors and dynamic performance can be 
determined by the simulation. However, it is not practical to lower the 


speed significantly because then the motor will be operating at a very low 


efficiency. 
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One area where additional work may be done is in utilizing the hybrid 
capabilities of the computer facility. Although the simulation has been 
done with only the Analog computer, the hybrid facility can be utilized to 
increase the effectiveness of the simulation. One very interesting study 
would be to perform an optimization study of the system. Even though 
prototype system performance is adequate, such a study would give a student 
the opportunity to learn the basics of hybrid computing and to study in the 
area of optimal control. 

The most important dynamic elements which were neglected were the fuel 
control valve and the voltage regulator. The governor control box stabilizer 
section was also neglected. Since the simulation required approximately 90. 
percent of the Analog hardware, it would be difficult to include reliable 
models of these components. A more detailed er of the generator set 
alone can be performed by replacing the induction motor load by a step load 
change. It will then be possible to include the effects of the fuel control 
valve, voltage regulator and stabilizing section. 

The generator set frequency regulation is extremely important for 
proper operation of the shipboard electronic equipment. As a result, the 
speed deviations of the gas turbine must be kept to a minimum when the 
generator is supplying power to the electronic equipment. A simulation of 
the generator set in response to step load changes can be used to design 
the governing system and to determine the governor settings required for 
best performance. 

The most direct benefit of the simulation would be to use it to demon- 
strate to operating 0 how the system performance can be altered by 


adjusting the governor settings. The desire to "leave well enough alone" 
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together with a lack of complete understanding of how a system operates 
will usually prevent operators from adjusting the 2 ۶۷۶۶ EO Cel better 
performance. The simulation can be used to instruct operators in how the 
various governor settings effect performance. With a better understanding 
of how the system behaves, it will be easier to adjust the governor set- 


tings for near optimum performance. 
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